可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

Orbit-attitude-vibration coupled dynamics of tethered solar power satellite

A new orbit-attitude-vibration coupled dynamic model of the tethered solar power satellite (Tethered SPS) is established based on absolute nodal coordinate formulation. The Hamilton’s equation of the system is derived by introducing generalized momentum through Legendre transformation. The correctness of the proposed model is verified by an example. The dynamic characteristics of the Tethered SPS are studied using the symplectic Runge-Kutta method. Simulation results show that the orbital radius and the total energy of the system are well preserved. The attitude of the system is unstable when the mass of the bus system is small. However, the attitude stability is dependent on some other parameters of the system, which requires further studies. It is also found that the average tether force/deformation can be roughly estimated by simplifying the solar panel as a particle. The proposed model can be used to study the orbit-attitude-vibration coupled dynamics and control problems.     

Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.     

逆金字塔式安全着陆点选择方法

针对现有地外天体软着陆安全着陆点选择方法的不足,设计了逆金字塔式安全着陆点选择方法,结合着陆器的落点精度和自主避障及地形坡度适应能力,建立了基于平整点、安全着陆点、可靠避障点和标称着陆点的地形安全性分析与定量搜索模型,该模型呈四层逆金字塔结构,每一层搜索的结果作为下一层搜索的输入,从弱约束到强约束逐步收敛至100%地形安全着陆区,嫦娥五号探测器飞行结果表明该方法能够快速确定安全着陆区及标称着陆点,并有效指导着陆下降过程的标称航迹设计和轨道控制策略制定。     

升力体式浮升混合飞艇多学科设计优化

升力体式混合飞艇是全球远距离大载重运输的重要选择，随着全球贸易的发展，逐渐成为国内外的研究热点。作为航空宇航技术、新能源技术和高性能材料技术相结合的新概念飞行器，混合飞艇设计过程需对多个学科进行综合考虑和优化。为了将多学科设计优化（MDO）方法引入到混合飞艇的总体设计中，将其分解为能源子系统、气动和推进子系统以及结构和重量子系统。在子系统模型构建的基础上，提出具有自适应能力的基于响应面的并行子空间优化（CSSO-RS）算法，将重量平衡和能量平衡作为实现远距离载重运输的约束条件，并提出爬升、日间巡航、滑翔和夜间巡航的多阶段任务剖面，以充分利用太阳能电池、燃料电池和锂电池的优势，实现混合飞艇的最优化设计。优化结果表明:具有自适应能力的优化算法在精确度和计算效率上均有明显的优势，同时重量分配的结果也为混合飞艇结构轻量化设计和能源系统设计提出了更高的要求。      

无人自转旋翼机建模与控制方法研究

以500kg级无人自转旋翼机为研究对象,研究了样例无人旋翼机的建模和自动飞行控制技术.首先,基于解析形式叶素积分法给出自转旋翼机的建模方法,并给出了样例旋翼无人机的数学模型;然后,设计了无人旋翼机的俯仰和滚转控制律,并在此基础上设计了速度、高度及侧偏控制律.最后,通过飞行试验验证了模型的准确性及控制律的有效性.     

中国空间探测领域40年发展

中国空间科学学会成立的40年，是中国空间探测逐渐走进世界舞台的40年，空间探测极大推动了空间科学和相邻学科的发展，也影响到经济、军事和日常生活诸多方面.本文简要回顾了从空间探测专业委员会成立的1980年至今，中国空间探测领域的主要发展历程，包括探空火箭、高空气球、科学卫星、月球与行星探测、载人航天空间探测、遥感卫星地面站等主要项目、进展和所取得的成果，对未来若干年空间探测的发展进行了展望.      

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

燃烧室头部激励的等离子体强化燃烧特性实验研究

等离子体助燃是一种新型的强化燃烧技术，近年来受到国内外学者的广泛关注。本文开创性地研制了基于旋转滑动弧等离子体的强化燃烧头部，建立了某型航空发动机三头部燃烧室实验件的等离子体助燃实验平台，验证了该等离子体强化燃烧技术应用于型号发动机燃烧室的可行性。实验研究等离子体助燃在不同余气系数和不同输入电压条件下对平均出口温度、燃烧效率、温度分布系数以及熄火边界的影响。实验结果表明，与正常燃烧相比，施加等离子体助燃后的燃烧效率有明显的提高，在输入电压为U0=240V，余气系数为 α=0.8时，等离子体助燃的燃烧效率提高3.24%。实施等离子体助燃后，燃烧室出口温度分布场分布得到明显的改善，在富油工况α=0.8，出口温度分布系数减少39.8%。等离子体助燃输入电压越高熄火边界扩展程度越明显，相比于正常工况条件下，等离子体助燃U0=240V的熄火边界扩宽了7.34%。     

Study of ionospheric D region changes during solar flares using MF radar measurements

During solar flares, the X-ray radiation suddenly increases, resulting in an increase in the electron density of the atmospheric D region and a strong absorption of short-wave radio waves. Based on Langfang medium frequency (MF) radar, this paper analyzed the variation characteristics of D region in the lower ionosphere from 62 km to 82 km. The analysis focused on multiple C-level and M-level solar flare events before and after the large-scale flare event at 11:53 (UT) on September 6, 2017. The results show that it is difficult to detect the electron density over 70 km in Langfang during solar flares, but the electron density value can be obtained as low as 62 km, and the stronger the flare intensity, the lower the detectable electron density height. Besides, the equal electron density height, the received power of X and O waves will also be significantly reduced during the flares, and the reduction of equal electron density height has a weak linear relationship with flare intensity.